Transverse momentum and rapidity dependence of HBT correlations 
in Au+Au collisions at y/s NN = 62.4 and 200 GeV 
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Two-particle correlations of identical charged pion pairs from Au+Au collisions at ^/s NN = 62 A 
and 200 GeV were measured by the PHOBOS experiment at RHIC. Data for the 15% most central 
events were analyzed with Bertsch-Pratt and Yano-Koonin-Podgoretskii parameterizations using 
pairs with rapidities of 0.4 < y T7r < 1.3 and transverse momenta 0.1 < kr < 1.4 GeV/c. The 
Bertsch-Pratt radii R and Re decrease as a function of pair transverse momentum, while R s is 
consistent with a weaker dependence. R and R 3 axe independent of collision energy, while Re 
shows a slight increase. The source rapidity y YKP scales roughly with the pair rapidity y TJr , indicating 
strong dynamical correlations. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Gz 



Recent experimental results from all four experiments 
at the Relativistic Heavy Ion Collider (RHIC) have con- 
cluded that Au+Au collisions at the top RHIC energy 
(%/ s jvn — 200 GeV) have produced an extremely hot and 
dense state of matter [jj, |2j, y, |4| . This matter may have 
degrees of freedom that are purely hadronic ("hadronic 
gas"), purely partonic ( "Quark-Gluon Plasma", QGP), 
or an admixture of both. Identical-particle correlation 
measurements (Hanbury-Brown and Twiss, HBT) yield 
valuable information on the size, shape, duration, and 
spatiotemporal evolution of the emission source. Because 
the dynamics of a hadron gas and a QGP are naively 
expected to be quite different, HBT may allow us to dis- 
criminate between these three scenarios pa, |(j . 

Experimentally, the correlation function C(q) is de- 
fined as 



C(q) 



P(pi,p 2 ) 
P( P i)P(p 2 ) 



(1) 



where pi and P2 are the particle four-momenta, 
P(pi,P2) is the probability of a pair being measured 
with relative four- momentum q = Pi — P2, and -P(pi) 
and -P(p 2 ) are the single-particle probabilities. The nu- 
merator is determined directly from data, while the de- 



nominator is constructed using a standard event-mixing 
technique. 

C(q) can be fit to the Bertsch-Pratt parameterization 
of a Gaussian source in three dimensions 0, 0, @ , 

C(q) = 1 + A e -(«o fl '+«' fl =+''< ii (+ 2 «»" ,i »<) (2) 

where qe is the component of q along the beam direction; 
q is the component along the pair transverse momentum 

kr = \{v T \ +PT2); an d <7s is the component orthogonal 
to the other two. The q qi cross-term vanishes only for 
symmetric collisions with acceptances centered around 
midrapidity. The A parameter represents the correlation 
strength and is expected to be unity for a completely in- 
coherent source. According to the definitions of q Q and q s , 
R Q probes a mixture of the spatial and temporal extent 
of the source, while R s measures only the spatial com- 
ponent. In the special case of a boost-invariant, trans- 
parent, azimuthally symmetric source, the ratio R Q /R S 
may be a good indicator of the duration of the emission 
of particles from the source. Predictions for this quantity 
from hydrodynamic and transport models varied by over 
an order of magnitude 0, El , but mostly focused on val- 
ues between 1.5-2.0, while the first results from RHIC at 
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y/s^ =130 GeV indicated a value near unity [lll Il2|. 
A more recent three-dimensional hydrodynamic calcula- 
tion (including opacity and tranverse flow) p| also was 
unable to reproduce the experimentally measured values 
for R a /R s . The results at y/s NN — 130 GeV also were 
consistent with a monotonic increase of Ri from AGS 
energies (y / s^" ~ 2-5 GeV), while R and R s remained 
roughly constant. It should be noted that it is possi- 
ble to fit the exp erimental data at RHIC to a variety 
of ansdtze [lj, [jjj, [lg , but the interpretation remains an 
open question. The wealth of HBT data over a large num- 
ber of variables and their experimentally-determined sys- 
tematic dependencies continue to create new challenges 
for theorists and add yet more pieces to the "HBT puz- 
zle". 

The data reported here for Au+Au collisions at 
t/s^T = 62.4 and 200 GeV were collected using the PHO- 
BOS two-arm magnetic spectrometer during RHIC Run 
IV (2004) and Run II (2001), respectively Details of 
the setup have been described previously |17| . The spec- 
trometer arms are each equipped with 16 layers of silicon 
sensors, providing charged particle reconstruction both 
outside and inside a 2 T magnetic field. The two sym- 
metric arms and frequent magnetic field polarity reversals 
allowed for a number of independent cross-checks in the 
two-particle correlation measurement. A two layer sili- 
con detector covering \rj\ < 0.9 and 25% of the azimuthal 
angle provided additional information on the position of 
the primary collision vertex. 

The primary event trigger was provided by two sets of 
16 scintillator paddle counters, which covered a pseudo- 
rapidity range 3 < \t)\ < 4.5. More information on PHO- 
BOS event selection and centrality determination can be 
found in [HHii. 

The details of the track reconstruction algorithm can 
be found in pol l2l| . Events with a reconstructed pri- 
mary vertex position between -12 cm < z vtx < 10 cm 
along the beam direction and -0.1 < x vtx < 0.2 cm and - 
0.05 < y v tx < 0.2 cm along the transverse directions were 
selected to optimize vertex-finding precision, track re- 
construction efficiency, and momentum resolution. Only 
particles which traversed the entire spectrometer were 
used in the analysis. A 3<r cut on the distance of clos- 
est approach of each reconstructed track with respect to 
the primary vertex {dca vtx < 0.35 cm) was used to reject 
background particles from decays and secondary inter- 
actions. The final track selection was based on the \ 2 
probability of a full track fit, taking into account mul- 
tiple scattering and energy loss. The momentum res- 
olution is Ap/p ~ 1% after all cuts. Particle identifi- 
cation was based on the truncated mean of the specific 
ionization loss in the silicon detectors. To identify pions, 
a cut of three RMS deviations away from the expected 
mean value of the specific ionization (dE/dx) for pions 
was applied. This is shown as the region between the 
solid black lines in Fig. 2] Possible contamination from 
other particle species was studied using HIJING 1.35 
and a GEANT 3.21 [23| simulation of the full detector, 
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FIG. 1: Distribution of the truncated average energy loss 
(dE/dx) as a function of reconstructed particle momentum. 
The region between the solid black lines corresponds to the 
±3<7 cut used to select pions for this analysis. 



applying the same reconstruction procedures to Monte 
Carlo events and data. The contamination from K + K + , 
K~K~, pp, pp, e + e + , and e~e - pairs is estimated to be 
less than 1%; non-identical pairs are estimated to con- 
tribute less than 10% throughout the entire fcy range. 
To reject ghost pairs, only one shared hit out of six in 
the weak-field region and two shared hits out of five in 
the strong-field region were allowed per pair. A two- 
particle acceptance cut was applied to both data and 
background; the criterion for pair acceptance was defined 
by A(j)+ 2 AO > 0.05 rad, where A(f> and AO are the rela- 
tive pair separation in azimuthal and polar angle, respec- 
tively. In this analysis, events in the most central 15% of 
the cross section were selected. HIJING was used to re- 
late the fraction of the cross section to (N part ), the mean 
number of participating nucleons [18( . The average num- 
ber of participants for these events was estimated to be 
(N part ) = 310 (303) at ^/s^ = 200 (62.4) GeV. About 
7.3 (3.3) million 7r + 7r + and 5.5 (3.0) million 7r~7r - pairs 
for y/s NN = 200 (62.4) GeV survive all cuts including 
centrality selection. 

The event-mixed background is constructed by com- 
bining single tracks from randomly selected different 
events into fake pairs. Events from widely spaced colli- 
sion vertices may yield a background that cannot repro- 
duce the same phase space as the actual data. This effect 
is particularly prominent for a small-acceptance detector 
such as PHOBOS. In this analysis, an "event class" is 
defined solely by the vertex resolution of the detector. 
Events were defined as belonging to the same event class 
if they had vertices within 0.025 cm along the beam axis 
and 0.05 cm along the vertical and horizontal axes. Back- 
ground events are only constructed from tracks mixed 
from the same event class. 

Systematic errors (90% C.L.) were determined by 
changing cuts in two-particle acceptance and single-track 
azimuthal acceptance, using different random seeds for 
mixed-event background generation, as well as varying 
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the event class definition to create background events 
from pairs with narrower and broader vertex ranges. Ad- 
ditionally, these studies were performed for each individ- 
ual spectrometer arm and the differences included in the 
systematic uncertainties. 

Because the background is constructed from tracks be- 
longing to different events, it does not a priori include 
multiparticle correlations (apart from a residual effect at 
the few percent level, which has been corrected following 
pil l25jl. In order to study the HBT correlation, it is 
necessary to apply a weight to account for the Coulomb 
effect. The Coulomb correction can be expressed solely 
as a function of the invariant relative 4-momentum (\i nv , 



F c (q mv ) _ /dr|^ c (r)| 2 g(r) 
F pl (q mv ) J df\iP pl {r)\ 2 S(f) 



(3) 



where S(r) is the distribution of the relative separation of 
the particle pairs at emission, Ri nv is the radius param- 
eter conjugate to Qi nv , and ip c and ip p i are the Coulomb 
and plane wave- functions, respectively. A closed- form 
approximation and numerical interpolation for this rela- 
tion was derived in [26l for A = 1. 
For variable A I27I. 



FRi„ v (linv, A) 



(l-A)+A(l + e- 



Fn im ,{<\ inv ) 



1 + Ae-* 



(4) 

This prescription was derived simultaneously and is 
nearly equivalent to the corrections applied by the 
CERES O, STAR HJ, and PHENIX experiments [H; 
our results showed no significant change using cither cor- 
rection method. The method is applied iteratively, suc- 
cessively fitting distributions of the correlation function 
C(c{inv) arid applying the fit values A and Ri nv to a new 
S(f). Typically 2 or 3 iterations are sufficient for conver- 
gence. 

The three-dimensional correlation functions were fit by 
Eq. (0) using MINUIT and the log-likelihood method. 
Table 0] shows the results of the fit for both tt + tt + and 
7r~7r~. The data were analyzed in the longitudinal co- 
moving system (LCMS) frame. In Fig. [3 the Bertsch- 
Pratt parameters and the ratio R Q / R s are presented as a 
function of kx for tt~t:~ pairs at ^s NN =62.4 GeV and 
200 GeV. The 200 GeV data are compared to data from 
STAR 29] and PHENIX H3- Not shown are values for 
A, which is roughly 0.5, and the cross-term R^, which is 
consistent with zero for all kx bins. The small vertical 
acceptance of the PHOBOS detector roughly translates 
into a small acceptance in q s , which is responsible for 
the large statistical and systematic uncertainties on the 
value of R s at low kx- Within errors, the values of R s 
are consistent with a weak dependence on kx, while the 
values of R and Ri decrease rapidly with increasing kx- 
Consequently, the ratio R /R s decreases from 1.24±0.19 
to 0.77 ± 0.11 for (k T ) from 0.2 to 0.7 GeV/c. It should 
be noted that the STAR and PHENIX data are mea- 
sured at midrapidity while the PHOBOS data are mea- 
sured at (y n7r ) = 0.9; additionally, STAR, PHOBOS, and 
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FIG. 2: Bertsch-Pratt parameters R , R s , Re and the ratio 
Ro/Rs for 7r~7r~ pairs at y?s NN = 200 GeV (left panel) and 
62.4 GeV (right panel) as a function of (fey). For comparison, 
data from STAR (open stars) and PHENIX [ajj (open 
crosses) are presented at 200 GeV. PHOBOS systematic er- 
rors (90% C.L.) are shown as boxes; systematic errors from 
STAR and PHENIX are not shown. 



PHENIX data are from the top 10%, 15%, and 30% of 
the cross-section respectively. 

Despite the very different experimental acceptances, 
the data at 200 GeV agree remarkably well. The trends 
a ^ -\/ s nn ~ 62.4 GeV are qualitatively and quantita- 
tively similar to those at 200 GeV, except that Re may 
be slightly smaller. 

The Bertsch-Pratt parameters A, R Q1 R S1 and Re are 
presented in Fig. as a function of center of mass en- 
ergy from ^/s NN = 2 to 200 GeV. The data are ir~ir~ 
pairs near midrapidity for comparable kx bins from nine 
different experiments. The data at ^/s NJV = 62.4 GeV 
begin to fill the large gap between the top SPS energy 
and the 130 GeV RHIC data. Although there is some dis- 
agreement among experiments, the data do not appear 
to exhibit any sharp discontinuities, but smoothly vary 
as a function of collision energy. The correlation function 
was also fit to the Yano-Koonin-Podgoretskii (YKP) pa- 
rameterization [31II32I 

C(q) = 1 + \ e -( q ^ B2 ^ 2 ( q u-P q ^ 2R l+^ 2 ^-P<ii0 2R2 ^ (5) 

where (3 is the longitudinal velocity of the source and 
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FIG. 3: Bertsch-Pratt parameters A, R 3 , R , and Rt as a 
function of ^/s NN for 7r _ 7r _ pairs. The presented data are 
near midrapidit y an d represent comparable kr bins from each 

experiment [ll 111 S ISl E3 lM El El E3 • PHOBOS 
data are represented by solid circles. Systematic errors are 
not shown. 



7 = l/yl — (3 2 , gj_ and gni are the relative 3-momentum 
difference projected in the transverse and longitudinal 
directions, respectively, and q T is the relative difference 
in energy. Employing this coordinate system provides the 
advantange of factorizing the velocity and duration of the 
source from the spatial parameters. In addition to fits in 
the LCMS frame, data were also fit in the laboratory 
frame as a cross-check and yielded consistent results. 

Table [H] shows the results of Yano-Koonin- 
Podgoretskii fits for ir + ir + and ir~ir~ at both energies. 
In Fig. the value of the source rapidity y YKP (as 
extracted from (3 in Eq. is plotted as a function of 
pair rapidity for n~n~ pairs with 0.1 < kr < 1.4 GeV/c. 
Open circles indicate the results reflected about midra- 
pidity. The data from NA49 H3 at ^/s^ = 17.2 GeV 
is also plotted; however, it should be noted the NA49 
data presented here cover only 0.1 < k T < 0.2 GeV/c. 
The source rapidity scales with the rapidity of the pair, 
indicating the presence of strong position-momentum 



correlations. A static source would exhibit no correlation 
and would correspond to a horizontal line (y YKP = 0). 
A source with strong dynamical correlations would 
correspond to a straight line along y YWP = y^. The data 
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FIG. 4: Dependence of the source rapidity y YKP vs. the pair 
rapidity y W7r for PHOBOS at 200 GeV (circles), PHOBOS at 
62.4 GeV (triangles), and NA49 (negatively charged hadrons) 
at 17.2 GeV |33| (squares). Open symbols represent data 
reflected about midrapidity. 



are consistent with the latter scenario: particles emitted 
at a given rapidity were produced by a source moving 
collectively at the same rapidity. In conclusion, we have 
measured HBT parameters in Au+Au collisions with 
the PHOBOS detector using both the Bertsch-Pratt and 
Yano-Koonin-Podgoretskii parameterizations at energies 
of = 62.4 GeV and 200 GeV. The ratio R /R s 

does not show significant deviation from unity over the 
entire fcr range. The data at 200 GeV show a nice agree- 
ment between three RHIC experiments with different 
acceptances, centrality measures, and rapidity, which 
implies that the source parameters are fairly insensitive 
to these variables at this energy. Additionally, the 
Bertsch-Pratt parameters evolve smoothly as a function 
of ,ys NN over nearly two orders of magnitude. 
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TABLE I: Bertsch-Pratt fit parameters for tt + 7v + and -it tv in the rapidity range 0.4 < y n „ < 1.3 for the 15% most central 
Au+Au collisions. Momenta are in GeV/c and radii are in fm. Errors are statistical (lcr) followed by systematic (90% C.L.) 

Species ^/s^ (k T ) A R a R s R t R 2 ol (fm 2 ) 

0.198 0.620±0.046±0.075 5.65 ± 0.33 ± 0.44 5.11 ± 0.69 ± 0.61 6.70±0.44±0.63 1.9 ±2.9 ±2.6 
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TABLE II: Yano-Koonin-Podgoretskii fit parameters for n + n + and f i in the kr range 0.1 < kr < 1.4 GeV/c for the 15% 

most central Au±Au collisions. Radii are in fm. Errors are statistical (lcr) followed by systematic (90% C.L.) 

Species ■ s fs^ (y^n) A R± R T (3 

0.602 0.669±0.057±0.074 5.05 ± 0.33 ± 0.30 5.75 ± 0.45 ± 0.70 0.0017±2.0± 1.4 -0.162±0.086±0.11 
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